home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Amiga Collections: Camelot
/
Camelot 098 (1990-12)(Swedish User Group of Amiga)(SE)(PD)[WB].zip
/
Camelot 098 (1990-12)(Swedish User Group of Amiga)(SE)(PD)[WB].adf
/
XLisp-Stat
/
Book
/
lik-gamma.lsp
< prev
next >
Wrap
Text File
|
1990-10-11
|
285b
|
14 lines
; book p.108
(require "functions/newton.lsp")
(def x (gamma-rand 30 4.5))
(def f (let ((s (sum (log x)))
(n (length x)))
#'(lambda (a)
(- (* (- a 1) s) (* n (log-gamma a))))))
(def df (make-derivative f .001))
(def ddf (make-derivative df .001))